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Abstract: We investigate the cosmological evolution in a universe governed by the 
extended, varying-mass, nonlinear massive gravity, in which the graviton mass is promoted 
to a scalar-field. We find that the dynamics always leads the varying graviton mass to 
q " zero at late times, offering a natural explanation for its hugely-constrained observed value. 

Despite the limit of the scenario towards standard quintessence, at early and intermediate 
times it gives rise to an effective dark energy sector of a dynamical nature, which can 
also lie in the phantom regime, from which it always exits naturally, escaping a Big- 
Rip. Interestingly enough, although the motivation of massive gravity is to obtain an IR 
modification, its varying-mass extension in cosmological frameworks leads rather to early 
and intermediate times modification, and thus to a UV modification instead. 
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1 Introduction 

The idea of adding mass to the graviton is quite old [1], but the necessary nonlinear terms 
[2] that can give rise to continuity of the observables [3, 4] lead also to Boulware-Deser (BD) 
ghosts [5], making the theory unstable. However, recently, a nonlinear extension of massive 
gravity has been constructed [6, 7] such that the Boulware-Deser ghost is systematically 
removed (see [8] for a review). The theoretical and phenomenological advantages, amongst 
which is the universe self-acceleration arising exactly from this IR gravity modification, 
brought this theory to a significant attention [9-49]. 

Despite the successes of massive gravity, in the case where the physical and the fiducial 
metrics have simple homogeneous and isotropic forms the theory proves to be unstable at 
the perturbation level [40], which led some authors to start constructing less symmetric 
models [13, 41]. However, in [50] a different approach was followed, namely to extend the 
theory in a way that the graviton mass is varying, and this was achieved by introducing an 
extra scalar field which coupling to the graviton potentials produces an effective, varying, 
graviton mass. 

In this work we desire to explore the cosmological implications of this "extended", 
varying-mass, massive gravity. As we show, at least in simple cosmological ansatzes, the 
dynamics always leads the varying graviton mass to zero at late times, and thus the theory 
has as a limit the standard quintessence paradigm. However, at intermediate times the 
varying graviton mass leads to very interesting behavior, with a dynamical effective dark 
energy sector which can easily lie in the phantom regime. Strictly speaking, although the 
motivation of massive gravity is to obtain an IR modification, its extension in cosmological 
frameworks leads rather to early and intermediate times modification, and thus to a radical 
UV modification instead. 
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2 Extended nonlinear massive gravity 



Let us briefly review the "mass- varying massive gravity" that was recently presented in [50] . 
Their construction is based on the promotion of the graviton mass to a scalar-field function 
(potential), with the additional insertion in the action of this scalar field's kinetic term and 
standard potential. Since such a modification is deeper than allowing for a varying mass, 
we prefer to call it "extended" nonlinear massive gravity. 
In such a construction the action writes as 



S = Id X\J —g 



M 2 _ „ „ „ , 1 



R + V(iP)(U 2 + a 3 U 3 + a A U A ) - -d^d^ - W(i/>) 



(2.1) 



where M p is the Planck mass, R the Ricci scalar, and ip is the new scalar field with 
W(ip) its standard potential and V{ip) its coupling potential which spontaneously breaks 
general covariance. Furthermore, as usual a 3 and «4 are dimensionless parameters, and 
the graviton potentials are given by 

u 2 = , u 3 = /cf^/c;, , u 4 = k\k? v w p K\ > ( 2 - 2 ) 

with Alj^/C^] = (fC^lC" — JCuJCp/2 and similarly for the other antisymmetric expressions, 
and 



= ^9^fABd p ^ A d v cp B . (2.3) 

As in standard massive gravity /ab is a fiducial metric, and the four 4> A (x) are the 
Stiickelberg scalars introduced to restore general covariance [51], and in the particular 
case where the Jab is the Minkowski metric they form Lorentz 4-vectors in the internal 
space and the theory presents a global Poincare symmetry, too. Finally, one can show that 
the above extended massive gravity is still free of the the Boulware-Deser ghost [50]. 

3 Cosmo-logical equations 

Let us now examine cosmological scenarios in a universe governed by the extended nonlinear 
massive gravity. For simplicity we consider the fiducial metric to be Minkowski 1 

Jab = VAB, (3.1) 

and without loss of generality we assume that the dynamical and fiducial metrics are diago- 
nalized simultaneously. For the physical metric we consider for simplicity a flat Friedmann- 
Robertson- Walker (FRW) form, although one can straightforwardly examine the non-flat 
case too [50]: 

d 2 s = -N(r) 2 dT 2 + airfSijdx^, (3.2) 

with N(t) the lapse function and a(r) the scale factor, and for simplicity for the Stiickelberg 
fields we choose the ansatz 

0° = b(r), ft = x\ (3.3) 



1 Note that this case includes the subclasses where Jab can be brought to the Minkowski metric by 
general coordinate transformation, as we can always choose a gauge for the Stiickelberg fields 4> A [50]. 
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Finally, in order to obtain a realistic cosmology one includes the usual matter action S m , 
coupled minimally to the dynamical metric, corresponding to energy density p m and pres- 
sure p m . Although the above specific application is only a simple subclass of the rich set 
of possible scenarios, it proves to exhibit very interesting cosmological behavior. 

Variation of the total action S+S m with respect to N and a provides the two Friedmann 
equations 



3Mj.iT 2 



PDE + Pr, 



-2MpH = p DE + PDE + Pm+ Pm, 



(3.4) 
(3.5) 



where we have defined the Hubble parameter H = a/ a, with a = da/ (Ndr), and in the end 
we set N = 1. In the above expressions we have defined the energy density and pressure 
of the effective dark energy sector as 



PDE 
PDE 



■tf + w(ip) + V(i/>) 



1 )[h(a) + fi(a)} 



1 



i> 2 - W(i>) - VW)h{a) - V^)bh{a) 



(3.6) 
(3.7) 



having also introduced the convenient functions 



Ma) 
Ma) 
Ma) 



3 h a 3 I 3 

a 



1 - - + a 3 1 

a 



3 - - + a 3 1 
a 



I) K)+-K)" 

a J 3 V a 



+ a 3 1 



(3.8) 



Note that there exists a - factor in the expressions of dark energy density and pressure, as 
well as in the functions /i, /2, f 3 and fa. The value 1 is the default reference that we have 
chosen in the fixed physical metric of form (3.2) and the Stiickelberg fields of form (3.3). 
One can set the reference value to an arbitrary one a re /, by rescaling a by a/a re f in the 
metric or transforming eft 1 to (p l /a re f in the Stiickelberg fields, but for simplicity throughout 
this work we set a re f = 1. 

One can easily verify that the dark energy density and pressure satisfy the usual 
evolution equation 



Pde + 3H(p D E + Pde) = 0, 

and we can also define the dark-energy equation-of-state parameter as usual as 

Pde 



w D e 



Pde 



(3.9) 



(3.10) 
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Note that in [50] the authors had named the aforementioned "dark energy" sector as 
"massive gravity" one, and the quantities pde and pde as pMG and pmg- However, since 
in this work we focus to late time cosmological behavior, we prefer the above name. 

Variation of the total action S+S m with respect to ip provides the scalar- field evolution 
equation: 



TT ■ dW dV 



--i][A(") + /i(«)] + 3&A(<0 

a 



0. (3.11) 

Furthermore, variation of S + S m with respect to b provides the constrain equation 

V^Hfria) + V(i>)f 2 (a) = 0. (3.12) 
Finally, one can also extract the matter evolution equation p m + 3H(p m + p m ) = 0. 

4 Cosmological behavior 

The cosmological implications of extended nonlinear massive gravity, namely equations 
(3.4), (3.5), (3.11) and (3.12), prove to be very interesting, however, at least in its present 
simple but general example, it can be radically different than the usual massive gravity. 
The reason for this is the constraint equation (3.12), which dramatically restricts V(ijj). 
Note that in standard massive gravity this equation imposes strong constraints on b(r) 
[13, 16, 49], but in the present model the extra freedom brings these constraints on V(V0 
instead, leaving b(r) free. In particular, for /2(a) 7^ (which is the case in general) we 
obtain 

V(4>(T)) = V e j ''to = ( a - l)[a 4 - (3q 3 + 2a 4 )a + (3 + 3a 3 + a 4 )a 2 ] ' ^ 

where we have used the definitions (3.8). The importance of the constraint equation (3.12) 
was not revealed in [50], where the authors claimed it to be a just a requirement for 
V(ip) 2 . However, from its solution (4.1) we can immediately observe that in general at 
late times V(iP(t)), that is the graviton mass, always goes to zero, independently of the 
details of the model such are the values of 03, «4, the potential W(ip), and the matter 
content. This means that the present scenario of extended nonlinear massive gravity, in 
cosmological frameworks, cannot provide the usual massive gravity 3 , and on the contrary 
it always gives the standard gravity along with the standard quintessence scenario [52, 53]. 
Similarly, once introduced, the scalar-field cannot be set to zero by hand, since this is not 
a solution of (3.11) and (3.12) (unless we also set V{ip) = but in this case the model 
coincides completely with standard quintessence), that is ip will always have a non-trivial 
dynamics. 

However, although at late times the present scenario coincides with standard quintessence, 
it can have a very interesting behavior at intermediate times. In particular, the dark energy 



2 In particular, the specific examples that were considered in [50] were exactly those fine-tuned parameter 
choices that lead to fi(a) = /2(a) = and thus to a trivial satisfaction of the constraint (3.12). 
3 In this sense perhaps the name "extended" nonlinear massive gravity is not appropriate. 
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sector is not only dynamical, but it can easily lie at the phantom regime [54-57] . This can 
be seen by observing pde and pde from (3. 6), (3. 7), which using the constraint equation 
(3.12) give 



Pde + Pde = ip - V(ip) b 



(4.2) 



So we can always find regions in the 03,04 parameter space, that can lead to pde + Pde < 
at some stage of the evolution (with a potential W(ifj) that will not lead to large ip), even 
if we require to always have pde > (which does not need to be the case in general). This 
null energy condition violation is always canceled at late times, where the vanishing of the 
graviton mass leads to wde > — 1- 

In order to present the above behavior in a more transparent way, we numerically 
evolve the scenario at hand choosing without loss of generality 6(r) = ln[a(r)] (with the 
dimensional coefficient set to 1 in Mp-units), for dust matter (w m = p m /p m = 0), using 
the redshift z = a^/a — 1 as the independent variable (with an the present scale factor), 
and imposing the present dark energy density Qde = PDE/(3MpH 2 ) to be « 0.72 and 
its initial value to be « 0. In Fig. 1 we present the effective dark-energy equation-of-state 
parameter woe as a function of z, in the case where the usual scalar field potential is 
W(tp) = Wo = const., and for two choices of the parameters 03, 04, Wo, Vq 4 . As described 
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Figure 1. Evolution of the dark energy equation-of-state parameter Wde as a function of the 
redshift z, in the case where the usual scalar field potential is W(ip) = Wq = const. The black-solid 
curve corresponds to 0:3 = —1, 0:4 = —1, Wq = 1, Vq = 40, while the red-dashed curve corresponds to 
CK3 = —2, «4 = —2, Wo = 2.5, Vo = 6, in Mp-units. The dotted -1-line is depicted for convenience. 



4 Note that the graviton mass is proportional to Vo but significantly downgraded by the denominator of 
(4.1), thus it is below the Planck mass even if Vo is chosen larger than Mp. 
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above, at early and intermediate times the coupling potential V(ip) is non-zero leading wde 
to exhibit a dynamical nature, which can lie in the quintessence regime (black-solid curve) 
or in the phantom regime (red-dashed curve). Additionally, as expected, since in both 
cases the usual scalar potential is constant, the system at late times, where the coupling 
potential V(i/j) of (4.1) becomes zero, tends to the cosmological constant behavior. 

In fig. 2 we present wde as a function of z, in the case where the usual scalar field 
potential has the exponential form W(ip) = Woe~ x ^, and for three choices of the parameters 
a3,at4,\,Wo,Vo. In all three sub-cases, wde presents a dynamical nature, which can be non- 
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Figure 2. Evolution of the dark energy equation- of- state parameter wde as a function of the redshift 
z, in the case where the usual scalar field potential takes the exponential form W(ip) = Wqb~ • 
The black-solid curve corresponds to = 2, = 2, A = 0.5, Wo = 15, Vq — 12, the red-dashed 
curve corresponds to a.3 = 1, (X4 = 1, A = 0.5, Wq = 5, Vo = 2.5, and the green-dotted curve 
corresponds to 03 = —2, 04 = —2, A = 1, Wq = 1, Vq = 2, in Mp-units. The dotted -1-line is 
depicted for convenience. In order to depict the limiting behavior towards quintessence, we have 
evolved the system towards the future and up to z = —0.9. 

monotonic or/and it can cross the phantom divide from below to above [58]. However, as 
described above, at late times, where the coupling V(ip) becomes zero, all three sub-cases 
are led to their usual quintessence limit, where the final wde is determined solely from the 
TV-potential exponent A [59], and such a limiting behavior can be clearly seen if we evolve 
the system towards the future and up to z = —0.9. 

In summary, as we can see the scenario at hand exhibits very interesting cosmological 
behavior at early and intermediate times, with a dynamical dark energy sector which 
can additionally lie in the phantom regime, before limit towards the standard quintessence 
scenario. Note that despite the phantom-crossing, at late times we always obtain wde > — 1 
since the vanishing of the graviton mass restores the null energy condition for the effective 
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dark energy sector, that is the universe will always escape from the phantom regime and 
the Big-Rip future [60, 61] that is common to the majority of phantom models. 

5 Discussion 

In this work we investigated the cosmological evolution in a universe governed by the 
extended, varying-mass, nonlinear massive gravity. Even for simple ansatzes the scenario 
proves to have a very interesting behavior, comparing with standard massive gravity. 

The first result is that the dynamics in cosmological frameworks always leads the 
varying graviton mass to zero at late times, and thus the theory possesses as a limit the 
standard quintessence paradigm. This is a great advantage of the present construction, 
since it offers a natural explanation of the tiny and hugely-constrained graviton mass 
that arises from current observations. The graviton mass does not have to be tuned to 
an amazingly small number, as it is the case in standard massive gravity, but it is the 
dynamics that lead it asymptotically to zero. 

Despite the vanishing of the graviton mass at late times, and the limit of the scenario 
towards standard quintessence, at early and intermediate ones it can lead to very interesting 
behavior. In particular, it can give rise to an effective dark energy sector of a dynamical 
nature, which can also lie in the phantom regime. The violation of the null energy condition 
for the effective dark energy sector at intermediate times arises naturally for suitable (not 
fine-tuned) regions in the Lagrangian parameters, and it is always canceled at late times due 
to the vanishing of the graviton mass. These features are in agreement with observations 
and they offer an explanation for the dynamical evolution of the dark-energy equation-of- 
state parameter, for its relaxation close or at the cosmological constant value, and also for 
the indicated possibility to have crossed the phantom divide. Moreover, even if it enters 
the phantom regime, the scenario at hand always returns naturally to the quintessence one, 
offering a solution to the Big-Rip fate of the the standard phantom scenarios. 

We mention here that although we performed the above analysis in a flat universe, 
and with the fiducial metric to be Minkowksi, qualitatively the obtained behavior is not 
a result of them (one can easily see that they hold in the non-flat case of [50] too) but it 
arises from the deeper structure of the theory, namely from the scalar-field coupling to the 
graviton potential. Thus, we do not expect the results to change in more general cases, 
unless one fine-tunes the theory. 

In conclusion, the extended, varying-mass, nonlinear massive gravity leads to very 
interesting cosmological behavior at early and intermediate times, while it always limits 
towards the standard quintessence scenario, where the graviton is massless and the extra 
scalar is only minimally coupled to gravity. Strictly speaking, although the motivation of 
massive gravity is to obtain an IR modification, its varying-mass extension in cosmological 
frameworks leads rather to early and intermediate times modification, and thus to a UV 
modification instead. 
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